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ABSTRACT 

The primary nature of the ^^C neutron source is very significant for the studies of the s-process 
nucleosynthesis. In this paper we present an attempt to fit the element abundances observed in 
16 s-rich stars using parametric model of the single neutron exposure. The calculated results 
indicate that almost all s-elements were made in a single neutron exposure for 9 sample stars. 
Although a large spread of neutron exposure is obtained, the maximum value of the neutron 
exposure will reach about 7.0 mbarn^^, which is close to the theoretical predictions by the 
AGB model. The calculated result is a significant evidence for the primary nature of the neu- 
tron source. Combining the result obtained in this work and the neutron exposure-initial mass 
relations, a large spread of neutron exposure can be explained by the different initial stellar 
mass and their time evolution. The possibility that the rotationally induced mixing process 
can lead to a spread of the neutron exposure in AGB stars is also existent. 
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1 INTRODUCTION 

The two neutron-capture processes, the slow (s-process) and the 
rapid (the r-process), occur under different physical conditions, 
and therefore are lik ely to arise in different astrophysical sites 
teurbidge et al.l [TQSTh . The most likely site for the s-process is 
the inter-shell region of a thermally pulsing AGB star, provided 
a suitable neutron source is active. Our present understanding of 
the behavior of s-process nucleosyn thesis has been reviewed by 
IBusso. Gallino & Wasserburi h999j) . There is a general consen- 
sus that the neutron source is the reaction ^'^C(q, n)^''0. In or- 
der to activate it, a partial mixing of protons (PMP) from the 
envelope down into the C-rich layers is required. PMP activates 
the chain of reactions (^^C(p, 7)"N(/3)"C(q, n)^''0). The phys- 
ical cause of PMP could be diffusive convective overshooting 
jHerwig et al] 1 19971; [ Goriely & Mowl avi 2000 ), rotationally in- 
duced mixing jLange r et al. 1999; Herwig, Lug aro & Langej2003l) 
and gravity waves iDenissenkov & Touli i2003h . The s-elements 
produced in the deep interior by successive neutron captures 
are subsequen tly broug ht to the surface by the third dredge- 
up ( iGallino et al] 1 199a) . Using the primary-like neutron source 
(^^C(a, n)^®0) and starting with a very low initial metallicity, 
most iron seeds are converted into ^'^*Pb. So, when third dredge up 
episodes mix the neutron capture products into the envelope, the 
star appears s-enhanced and lead-rich. Therefore, if the standard 
PMP scenario holds, all s-process-enriched AGB stars with metal- 
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licities [Fe/H]sC-1.3 are thus predicted to be Pb stars ([Pb/hs]jsl, 
where "hs" denotes the 'heavy' s-process ele ments such as Ba, La, 
Ce), independent of their mass and metallicity jGorielv & Mowlavil 
200(|. 

The first three such lead st ars (HD187861, Hp 224959, 
HD 196944) have been repor ted by IVan Eck et al.1 ( 1200 ll) . At the 
same time. lAoki et al. I ( l200lh found that the slightly more metal- 
deficient stars LP 625-44 and LP 706-7 are enriched in s-elements, 
but cannot be considered as lead stars, in disagreement with the 
standard PMP predictions. 

The large observation data spreads of [Pb/hs] are strong in- 
dication to suspect a large intrinsic spread of integrated neutron 
irr adiations. A large sp read of ^"^C pocket efficiencies is proposed 
bv lStraniero et alj j2005» in order to explain the spreads of [Pb/hs]. 
It should be stressed here that the primary nature of the ^ ^C neu- 
tron source are rather robust toorielv & Mowlavil |2000|) . In the 
framework of the PMP scenario, there is no obvious degree of 
freedom that could be used to r educe the ^^C pocke t efficien- 
cies in low-metallicity AGB stars jVan Eck et id]|2003 '). The pos- 
sibility that the rotationally induced mixing process may lead to 
a spread of the neutron exposu re in AGB stars is proposed by 
iHerwig. Lugaro & Langej j2003l) . lsiess. Gorielv & LangeJ l [20oi) 
explored the effects of rotationally induced mixing on the nucle- 
osynthesis of s-elements during the TP-AGB phase. The results 
show that rotation of the AGB star quenches the s-process effi- 
ciency because of the contamination of the ^''C layer by ^*N. They 
find that although rotational mixing is an efficient mechanism to 
trigger the third dredge-up, but this process fails to reproduce the 
observed strong s-process overabundances. The dependence of the 
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induced mixing on the initial value of the rotational velocity is also 
largely unknown. So the fundamental problems, such as the for- 
mation and the consistency of the ^^C pocket, the development of 
the third dredge-up and the neutron exposure signature in the in- 
tetpulse, current ly exist in the models of AGB stars. In contrast 
to other studies, ICui & Zhand ( |2006|) find that a large spread of 
the ^^C efficiency is not needed to explain the observed spread of 
[Pb/hs], but this comes naturally from the range of different ini- 
tial st ellar masses and their time evolution (see also Bonacic et al., 
l2006h . The neutron exposure per circle deduced for the s-rich star s 
lies between 0.45 and 0.88 mbarn"^ jZhang. Ma & Zhoull200d) . 
For the metallicites of the Pb-enhan ced stars, based on the primary 
natur e of the ^^C neutron source dBusso. Gallino & Wasserbura 
Il999i) . th e neutron exposure per interpulse will reach about 7.0 
mbam" dCui & Zhan j|2006l) which are about 10 times of the re- 
sults obtained by Zhang. Ma & Zhoul ( l2006h for s-rich metal-poor 
stars. 

Another possibility is that the s-process material has ex- 
perienced only a few neutron exposures in the convective He- 
burning shell. This is consistent with a proposed mechanism 
for the s-process in metal-poor AGB stars with [Fe/H]<-2.5 
dpuji moto. Ikeda & Iben 2000). These authors proposed a scenario 
in which the convective shell triggered by the thermal runaway de- 
velops inside the helium layer. Once this occurs, ^^C captures pro- 
ton to synthesize ^^C and other neutron- source nuclei. The ther- 
mal runaway continues to heat material in the thermal pulse so 
that neutrons produced by the ^^Ne(a, n)^^Mg reaction as well 
as the ^^C(a, n)^^0 reaction may contribute. It is possible that 
only one episode of pr oton mixing into He intershe l l layer oc- 
curs in metal-poor stars jpuiimoto. Ikeda & Ibenll200(il : lAoki et aP 
l200ltllwamoto et al. '2002). After the first two pulses no more pro- 
ton mixing occurs although the third dredge-up events continue 
to repeat, so the abundances of the s-rich metal-poor stars can be 
characterized by only one neutron exposure. Detailed stellar evo- 
lution calculations are therefore highly desired, in order to clarify 
which site is the most likel y to dominate the s -process in metal- 

i ioor AGB stars [interpulse i Gallino et al.lll99^ . or thermal pulse 
Fuiimoto, Ikeda & Ibenll2000h l. 

Obviously, the detailed study of s-rich stars are needed in 
order to make progress in our understanding of the s-rich phe- 
nomenon, investigate what its physical reasons might be and for 
constraining what the possible physical conditions are. These rea- 
sons motivated us to start a systematic study for s-rich stars. In this 
paper we use the parametric approach of only one neutron expo- 
sure to investigate the characteristics of the nucleosynthesis path- 
way that produces the abundance ratios of s-rich objects. 



2 PARAMETRIC MODEL OF THE SINGLE NEUTRON 
EXPOSURE 

In order to investigate the efficiency of the s-process, the elemen- 
tal abundances of s-rich stars are particularly useful. There have 
been many theoretical studies of s-process nucleosynthesis in low- 
mass AGB stars. Unfortunately, however, the precise mechanism 
for chemical mixing of protons from the hydrogen-rich envelop 
into the ^^C-rich layer to form ^^C-pocket (i2C(p,7)^^N(e+,i/)^^C 
) is still unknown. In this paper we analyze the direct observa- 
tional constraints provided by the photospheric composition of s- 
rich stars. Preliminary attempts to fitting their abundances eithe r 
using parameterized s-process distributions (e.g. lAoki et al.ll200ll ) 
give encouraging results. For this reason, we use the parametric ap- 



proach of only one neutron exp osure, with many of neutron-capture 
rates updated teao et al.ll200a) , to investigate what physical con- 
ditions are possible to reproduce the observed abundance pattern 
found in the s-ric h stars tAoki et al. 200 1, 20()2; Cohen et al 200| ; 
Johnso n & Boltd| 2002l; ULuc atello et al.l l2003l ; iBarbuv et al] 12003 : 
Ivans et al..2005l : 



Jonsell et al..,2006') . 



We explored the origin of the neutron-capture elements in s- 
rich stars by comparing the observed abundances with predicted s- 
and r-process contribution. The ;-th element a bundance in the en- 
velop e of the star can be calculated as follows jZhang. Ma & Zhoj 
l2006l) : 



Ni{z) = c»iv,, , + ,.10 



[Fe/H] 



(1) 



where Z is the metallicity of the star, Ni, s is the abundance of the 
;-th element produced by the s-process in AGB star and Ni^ the 
abundance of the ;-th element produced by the r-process(per Si=10^ 
at Z=Z0), Ca and Cr are the component coefficients that corre- 
spond to relative contribution from the s-process and the r-process. 
It should be noted that the s-process abundance in the envelope of 
the stars could be expected to be lower than the abundance pro- 
duced by the s-process in AGB star because the material is mixed 
with the envelopes of the primary (former AGB star) and secondary 
stars. 

For the single neutron exposure model, the overlap factor r, 
which is the fraction of material that remains to experience sub- 
sequent neutron exposures, is not a parameter because there is no 
material that experience subsequent neutron exposures. So there are 
only three parameters in the parametric model. They are the neu- 
tron exposure Ar, the component coefficient of the s-process Cs 
and the component coefficient of the r-process Cr. 

We can carry out s-process nucleosynthesis calculation of sin- 
gle neutron exposure to fit the abundance profile observed in the 
s-rich stars, in order to look for the minimum in the three- 
parameter space formed by Ar, Cs and C,.. Using the method pre- 
sented by Aoki et al. (2001), we chose Sr as the representative for 
the first peak elements, Ba as the representative for the second peak 
elements and Pb as the representative for the third peak elements, so 
the uncertainties of the parameters are determined by the error lim- 
its of the representative elements. The adopted initial abundances of 
seed nuclei lighter than the iron peak elements were taken to be the 
solar-system abundances, scaled to the value of [Fe/H] of the star. 
Because the neutron-capture-element component of the interstellar 
gas that formed very mental-deficient stars is expected to consist 
of mostly pure r-process elements, for the other he avier nuclei we 
use th e r-process abundances of the solar system jArlandini et al.l 
1 19991) . normalized to the value of [Fe/H]. The adopted abundances 
of r-process nuclei in equation ( 1 ) are taken to be the solar-system r- 
process abundances ( Arlandini et al. 1999) for the elements heavier 
than Ba, for the other lighter nuclei we use so lar-system r-process 
abundances multiplied by a factor of 0.4266 jZhang. Ma & ZhotJ 
l2006l) . 



3 RESULTS AND DISCUSSION 

With the observed data in 16 sample stars, the model parameters 
can be obtained from the parametric approach. The results of the 
neutron exposures Ar, Cs, Cr and s-process fractional contribu- 
tions for Ba and Eu, fsa.s and fEu.s are listed in the table 1. For 
9 sample stars (the first 9 stars listed in table 1), a single neutron 
exposure fits well the observed data within the error limits of the 
representative elements, whereas for others, although we can find 
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the minimum in the three-parameter space, a single neutron ex- 
posure does not provide the parameters within the error limits for 
the entire representative elements. 

Figure 1 shows our calculated best-fit results for the 9 sam- 
ple stars. In order to facilitate the comparisons of the theoretical 
abundances with observations, the observed abundances of heavy 
elements are marked by filled circles in these figures. We note from 
the figure that, for most stars, the curves produced by the model are 
consistent with the observed abundances within the error limits. 
The agreement of the model results with the observations provides 
a strong support to the validity of the parametric model adopted in 
this work. 

The neutron e xposure deduced fo r s-rich stars lies between 0.7 
and 7.12 mbarn~^ . lAoki et^ ( 120011) have reported a neutron ex- 
posure, At ~ 0.71mbarn~^ for metal-poor star LP 625-44 and 
At - O.Smbarn,-^ for LP 706-7, our calculated results are close 
to their values. iGallino et al. 

I |l998) pointed out that the neutron 
density is relatively low, reaching lO'^cm"'' at solar metallicity, cor- 
responding to At ~ 0.2mbarn~^ . Since the ^^C neutron source 
is of primary nature, the typical neutron density in the nucleosyn- 
thesis zone scales roughly as l/Z"", from Zq down to 1/50 Zq. 
At lower metallicities, the effect of the primary poisons prevails 
( IBusso, Gallino & Wasserbu^ll999l) . For the metallicities of the 
sample stars([Fe/H]<-2.0), based on the primary nature of the ^''C 
neutron source, the mean value of neutron density in t he ^'^C-pocket 
is abo ut 35 times of the value for solar metallicity (Gallino et al. 
1 19991) . which indicates that the neutron exposure per interpulse will 
reach about 7.0 mbarn^^. It is interesting to note that, although a 
large spread of neutron exposure is shown in table 1, the maximum 
neutron exposure is close to this value. 

We discuss the uncertaint y of t he parameters using the 
method presented by lAoki et al] ( 1200 ih . Figure 2 shows the cal- 
culated abundances loge(Pb), loge(Ba) and loge(Sr) as a func- 
tion of the neutron exposure At in a model with Cr=10.7 and 
Cs=0.0014. These are compared with the observed abundances of 
HE 0024-2523. There is only one region in Figure 2, At=7.12'^q\q 
mbarn"^, in which all the observed ratios of four representative 
elements can be accounted for within the error limits. The bot- 
tom panel in Figure 2 displays the reduced value calculated 
in our model with all detected elemental abundances being taken 
into account and there is a minimum, with x^=4.27, at At = 
7.12mbarn~^ . For most sample stars shown in Figure 1, the uncer- 
tainties of the neutron exposure are smaller than that of HE 0024- 
2523. 

It is worth commenting on the behavior of log£(Sr), loge(Ba) 
and loge(Pb) as a function of the neutron exposure At seen 
in Figure 2. The nonlinear trends displayed in the plot reveal 
the complex dependence on the neutron exposure. The trends 
can be illustrated as follows. Starting at low neutron exposure 
and moving toward higher neutron exposure values, they show 
how the Sr peak elements are preferentially produced at nearly 
At OAmbarn'^ . At larger neutron exposure (e.g.. At ~ 
0.7mbarn~^), the Ba-peak elements become dominant. Then a 
higher value of loge{Pb) ^ 2 follows at At — 1.5mbarn^^ . 
In this case, the s-process flow extends beyond the Sr-peak and 
Ba-peak nuclei to cause an accumulation at ^"**Pb. At very high 
neutron exposure values (At>2 mbarn^^), the role of Fe nuclei 
as seeds of the s-process is replaced partly by lighter (intermedi- 
ate atomic mass) nuclei. Neutron capture on them can cross the 
iron peak, thus allowing the s-processing on heavy isotopes to con- 
tinue. It is clear from Figure 2 that there is another maximum, with 
loge{Ba) — 1.0 at At — 5.5mbarn~^ . In this case, the Ba-peak 



elements become dominant again. For extremely high neutron ex- 
posure values (At>6 mbam~^), the lighter seeds are also con- 
verted into ^"**Pb. Clearly, loge(Pb) is very sensitive to the neutron 
exposure. 

Our model is based on the observed abundances of the s-rich 
stars and the nucleosynthesis calculations, so the uncertainties of 
those observations and measurement of the neutron-capture cross 
sections will be involved in the model calculations. For HE 0024- 
2523, abundances of four neutron capture elements (Sr, Ba, La and 
Pb) and upper limits f or an additional three elements (Y, Zr and Eu; 
iLucatello et"Zll2003l) have been taken into account in our model. 
The three upper limits will enlarge the uncertainties of the param- 
eters. We note from Table 1 that for four stars (CS 31062-50, HE 
2148-1247, LP 625-44, HE 0024-2523), the reduced are larger 
than 2. The probability that x^ could be this large as a result of 
random errors in the measurement of the neutron-capture cross sec- 
tions and abundances of the neutron-capture elements is less than 
2%. We find that all these uncertainties can not explain the larger 
errors of neutron-capture elements, such as Zr in HE 2148-1247 
and Y in LP 625-44. This implies that our understanding of the true 
nature of s-process or r-process is incomplete for at least some of 
these elements (Travaglio et al. 2004). 

It was possible to isolate the contributions corresponding to 
the s- and r-process by our parametric model. In the Sun, the ele- 
mental abundances of Ba and Eu consist of significantly different 
combinations of s- and r-process isotope contributions, wi th s:r ra- 
tios f or Ba and Eu of 81:19 and 6:94, respectively t Arlandini et al.l 
Il999l) . The Ba and Eu abundances are most useful for unraveling 
the sites and nuclear parameters associated with the s- and r-process 
correspond to those in extremely metal-poor stars, polluted by ma- 
terial with a few times of nucleosynthetic processing. We explored 
the contributions of s- and r-process for these two elements in the 
s-rich stars. In table 1 we display the s-process fractions calculated 
from equation (1) for the 16 sample stars. Clearly, for the first 9 
stars listed in table 1, the s:r ratios for Ba and Eu are larger than 
90:10 and 15:85, which are larger than the ratios in the solar sys- 
tem. The abundances of r-elements, such as Eu, in s-rich stars are 
usually higher than those in normal stars, so s-rich stars seemed 
to be also enriched in r-elements, although in a lower degree than 
s-elements. 



4 CONCLUSIONS 

In this work, s- and r-process in s-rich stars were studied using 
the parametric approach of single neutron exposure. Theoretical 
predictions for abundances starting with Sr fit well the observed 
data for 9 sample stars (the first nine stars listed in table 1), pro- 
viding an estimation for neutron exposure occurred in AGB stars. 
The calculated results indicated that almost all s-elements were 
made in the first neutron exposure for 9 sample stars. Once this 
happens, after only one time dredge-up, the observed abundance 
profile of the s-rich stars may be reproduced in a single neu- 
tron exposure. This should be consistent with a proposed mecha- 
nism for the s-process in meta l-poor AGB stars with [Fe/H]<-2.5 
dFuiimoto, Ikeda & Ibenlbood) . 

For the third dredge-up and the AGB model, several impor- 
tant properties dep end primarily on the core mass (Iben 197^ 
iKarakas et ai]|2002h . The analytical formula for AGB stars given 
by iibd iiaTi" sh ow that the overlap facto r decreases with in- 
creas ing core mass. lKarakas et al. I ( |2002|) and lHerwig et all ( l2OO0l 
|2004 have also found that the third dredge-up is more efficient 
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Table 1. Observed abundance ratios and the derived parameters for s-rich stars 



Star 


[Fe/H] 


[Ba/Fe] 


[Eu/Fe] 


[Pb/Ba] 


Ar 
(mbarn^^) 


Cs 


Cr 


fBa,s 

% 


fSu.s 
% 




CS 29526-110 


-2.38 


2.11 


1.73 


1.19 


4.10 


0.0050 


49.6 


89.5 


14.5 


0.726588 


CS 22898-027 


-2.25 


2.23 


1.88 


0.61 


0.88 


0.0035 


62.8 


90.6 


23.3 


1 .435540 


CS 31062-012 


-2.55 


1.98 


1.62 


0.42 


0.84 


0.0017 


34.7 


91.1 


22.9 


1.162923 


CS 31062-050 


-2.32 


2.30 


1.84 


0.60 


0.82 


0.0041 


49.3 


94.9 


34.0 


2.229041 


HE 2148-1247 


-2.30 


2.36 


1.98 


0.76 


0.90 


0.0049 


64.5 


92.3 


28.1 


2.130968 


LP 625-44 


-2.71 


2.74 


1.97 


-0.19 


0.70 


0.0058 


63.6 


96.7 


35.4 


3.006556 


LP 706-7 


-2.74 


2.01 


1.40 


0.27 


0.84 


0.0019 


15.9 


96.1 


41.3 


0.602440 


HE 0024-2523 


-2.71 


1.46 


1.10 


1.84 


7.12 


0.0014 


10.7 


92.2 


19.7 


4.273823 


HE 0338-3945 


-2.42 


2.41 


1.94 


0.69 


0.90 


0.0055 


55.5 


94.0 


33.8 


1.435311 


CS 22948-27 


-2.47 


2.26 


1.88 


0.61 


0.76 


0.0026 


57.6 


92.7 


22.9 


0.990173 


CS 29497-34 


-2.90 


2.03 


1.80 


0.92 


3.80 


0.0035 


56.9 


82.6 


8.5 


2.399216 


HD 196944 


-2.25 


1.10 


0.17 


0.80 


3.30 


0.0010 


0.5 


99.2 


69.0 


2.598222 


CS 30301-015 


-2.64 


1.45 


0.2 


0.25 


3.76 


0.0007 


0.7 


98.6 


58.3 


3.644140 


CS 29497-030 


-2.57 


2.32 


1.99 


1.33 


4.54 


0.0053 


86.3 


85.2 


10.3 


2.459433 


CS 30322-023 


-3.41 


0.54 


-0.51 


0.95 


3.92 


0.0002 


0.1 


99.3 


72.2 


1.576791 


CS 22183-015 


-3.12 


2.09 


1.39 


1.08 


1.28 


0.0056 


21.5 


85.0 


17.0 


1.291154 




40 50 60 70 80 40 50 60 70 80 40 50 60 70 80 

Z {Proton Number) Z (Proton Number) Z (Proton Number) 



Figure 1. The Best fit to observational results of metal-deficient stars. The black circles with appropriate error bars denote the observed element abundances, 
the solid lines represent predictions from s-process calculations considering r-process contribution. 



for the AGB stars with larger core masses. Since the core mass of 
AGB star at lo w metallicity is larger than high metallicity obviously 
( IZiilstral2004l) . the overlap factor should be significantly smaller at 
low metalli city. This is consis t with the small overlap factor, r=0. 1, 
deduced bv lAokietal.) bOOlh for metal-poor stars LP 625-44 and 



LP 706-7. In an evolution model of AGB stars, a small r may be also 
realized if the third dredge-up is deep enough for s-processed ma- 
terial to be diluted by extensive admixture of unprocessed material. 
Clearly, the behavior of very deep dredge-up (r=0) is also single 
neutron exposure. Once this happens, no matter how many pulses 
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Figure 2. Calculated abundances loge(Pb), loge(Ba),log£(Sr) and reduced as a function of the neutron exposure per pulse, At, in a model with Cr=10.7 
and Cs=0.0014. Solid curves refer to the theoretical results, and dashed horizontal lines refer to the observational results with errors expressed by dotted lines 
for HE 0024-2523. The shaded area illustrates the allowed region for the theoretical model. 



may follow, the observed abundance profile of the s-rich stars may 
be reproduced in several similar single neutron exposures. 

For the sample stars, the maximum value of the neu- 
tron exposure obtained in this work reach about 7.0 mbarn"^, 
which is close to the theoretical pr edictions by the AGB model 
( IBusso, Gallino & Wasserbu^ll999f) and the maximum value of 
the case of r adiative ^^C-burning in the low mass AGB stars 
jCui & Zhan j |2006). So the calculated result should be a signifi- 
cant evidence for the primary nature of the neutron source. 

The large spread of neutron exposure given in table 1 could 
be explained when varying the initial mass that affect dredge- 
up events and nucleosynthetic processes along the AGB evolution 
dCui & Zhansll2006l : lBonacic et al.l2006l) . Certainly, the possibility 
that the rotationally induced mixing process may lead to a spread of 
the neutron exposure in AGB stars is also existent. For intermediate 
mass stars, the rapid rotation and significant amount of shear at the 
bottom of the convective envelope in AGB stars are hard to avoid. 
The rotations of the AGB star will quench the s-process efficiency 
and reduce the neutron exposure. Through a distribution of initial 
rotation rates, it may be lead to a natural spread of neutron exposure 
in AGB stars. In table 1, one interesting result is that, based on the 
values of the neutron exposure, the sample stars should be divided 
into two groups: one group (include 8 sample stars) concentrate on 
the value of about 0.8 mbam ~^ which correspond to the value of 
2.5M0 AGB star calculated bv lCui & Zhand ( |2006|) , and the others 



(include 8 also sample stars) are in the range of 1.3—7.12 mbarn"^ 
which correspond to the range of 1.5—2.0 Mq. Should the observed 
abundance pattern of the group one belongs to be polluted by the 
rapid rotators or the AGB stars with initial mass of about 2.5—3.0 
Mq? Obviously, large uncertainties still remain in this topic. More 
in-depth theoretical and observational studies of s-rich stars will 
reveal the characteristics of the s-process at low metallicity, such 
as their initial rotation rate and initial mass dependence, and the 
history of enrichment of s- and r-elements in the early Galaxy. 
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